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ABSTRACT 


This  study  develops  a  method  for  analytically  propagating  a  covariance  matrix 
to  a  maneuver  condition  to  be  used  in  linear  covariance  analysis  for  planning  the 
rendezvous  phase  of  a  space  mission.  With  the  generalized  formulation  of  a 
condition  transition  matrbc,  an  analytic  method  of  propagating  an  augmented 
covariance  matrix  to  any  scalar  terminal  maneuver  condition  is  presented.  The 
twenty-sbc  dimensional  augmented  covariance  matrbc  used  in  this  study  includes 
navigation  state  errors,  state  dispersions,  and  time  errors  for  both  the  chaser  and 
target  craft. 

The  method  is,  first,  analytically  developed.  The  vehicles  are  brought  to  the 
desired  rendezvous  condition  by  linearizing  the  motion  at  the  maneuver  condition 
point  and  allowing  the  time  of  flight  to  vary  slightly.  The  analytic  propagation 
technique  is  then  validated  by  comparison  to  a  stochastic  Monte  Carlo  simulation 
for  the  case  of  several  elevation  angle  conditions  which  might  be  used  to  trigger  an 
initial  rendezvous  intercept  burn.  The  validity  of  linearizing  the  motion  about  the 
terminal  point  is  substantiated  with  the  same  simulation. 
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CHAPTER  1 
INTRODUCTION 


TTie  purpose  of  (his  study  is  to  develop  an  analytic  method  for  propagating 
a  covariance  matrix  to  a  maneuver  condition  to  be  used  in  linear  covariance  analysis 
for  planning  the  rendezvous  phase  of  a  space  mission.  Rendezvous  performance 
analysis  usually  includes  the  use  of  linear  covariance  and  deterministic  simulation 
techniques.  Most  rendezvous  profiles  include  maneuvers  to  conditions  (eg.:  the 
position  of  the  chaser  craft  relative  to  the  target),  as  well  as  the  more  standard 
maneuvers  to  a  fixed  time.  This  study  presents  an  analytic  method  of  propagating 
the  covariance  matrix  to  a  condition  other  than  time. 

This  chapter  motivates  and  introduces  the  subject:  Section  1.1  provides  the 
motivation  for  developing  the  method,  Section  1.2  details  a  brief  background  of  the 
current  approach  to  propagating  to  an  elevation  angle  condition.  Section  1.3 
discusses  previous  work  that  forms  a  starting  point  for  this  study,  and  Section  1.4 
overviews  the  remaining  chapters. 
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1.1  Motivation 


Future  missions  to  Mars  will  stretch  current  theories  and  practices  in  mission 
planning  and  operations  well  beyond  those  used  to  go  to  the  moon.  Several  factors 
serve  to  greatly  complicate  the  problem.  Probably  the  greatest  complicating  factor 
is  the  distance  involved.  Because  of  the  distance,  most  operations  in  the  vicinity  of 
Mars  will  need  to  be  autonomous.  Recent  mission  architectures,  whether  manned 
or  robotic,  have  pointed  to  a  need  for  an  elliptic  rendezvous  capability.  Great  fuel 
savings  can  be  gained  by  leaving  the  heavy  interplanetary  ship  in  a  high  energy 
elliptic  orbit.  The  cost  of  taking  only  the  landing  craft  down  to  a  lower  orbit  and/or 
the  surface,  and  then  lifting  it  back  to  the  high  energy  orbit  would  be  much  smaller 
than  the  cost  of  putting  the  interplanetary  ship  in  a  low  circular  orbit.  This  mission 
scenario,  coupled  with  the  remoteness  of  Mars,  leads  to  a  more  complicated 
rendezvous  problem  than  the  near  circular  rendezvous  typical  of  Apollo  and  current 
Space  Shuttle  operations. 

Some  of  the  questions  that  need  to  be  re-addressed  in  an  elliptic  rendezvous 
problem  are  [3]: 

•  the  shape  and  size  of  the  catch-up  orbit, 

•  where  in  the  orbit  the  rendezvous  should  be  initiated, 

•  the  condition  on  which  to  trigger  the  intercept  burn, 

•  where  in  the  orbit  the  rendezvous  should  occur, 

•  the  geometry  of  the  final  approach  to  the  target,  and 

•  whether  the  rendezvous  scheme  can  be  made  tolerant  of  navigation 
uncertainties. 
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All  these  questions  are  interrelated,  and,  clearly,  some  trade-offs  will  have  to  be 
made.  Near  circular  rendezvous  experience  provides  a  starting  point,  but 
complications,  like  the  orbital  velocity  vector  no  longer  being  aligned,  essentially,  in 
the  horizontal  direction,  may  dictate  different  solutions  than  those  currently  used  for 
near  circular  orbits. 

Because  of  the  trade-offs  involved,  mission  planning  is  an  iterative  process. 
The  analytic  method  set  forth  in  this  study  was  devised  to  simplify  and  increase  the 
accuracy  of  the  search  for  answers  to  these  questions  during  the  mission  planning 
phase  by  providing  an  analytic  method  of  propagating  the  covariance  matrix  to  a 
condition.  While  the  method  of  propagating  to  a  condition  will  be  completely 
general,  this  study  will  focus  on  observable  conditions  which  might  be  used  as  the 
trigger  point  for  the  intercept  burn. 


1.2  Background 

Presently,  the  condition  used  to  trigger  the  intercept  burn  for  co-circular 
rendezvous  profiles  and  for  reducing  approach  trajectory  dispersions  for  shuttle 
rendezvous  operations  is  the  elevation  angle  of  the  line  of  sight  from  the  chaser  to 
the  target  craft,  measured  from  the  chaser’s  local  horizontal.  If  the  elevation  angle 
is  held  constant,  then  for  all  dispersed  trajectories,  state  errors  perpendicular  to  the 
line  of  sight  between  the  spacecraft  will  be  zero  at  the  condition.  The  current 
method  of  propagating  a  covariance  matrix  to  this  condition  uses  an  empirical 
approach.  The  matrix  is  first  extrapolated  to  the  time  at  which  the  nominal  orbit 
meets  the  condition.  It  is,  then,  empirically  ’shaped’  so  the  covariance  reflects  the 
zero  error  condition  perpendicular  to  the  line  of  sight  and  the  resulting  dispersions 
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that  exist  at  the  condition  as  seen  from  deterministic  simulations.  Monte  Carlo 
simulation  results  confirm  this  method  is  essentially  equivalent  to  extrapolating  to 
the  condition.  This  study  seeks  to  find  an  analytic  method  of  extrapolation  to  the 
condition  by  allowing  for  the  possibility  of  time  slips  which  are  necessary  to  meet  the 
condition  in  the  presence  of  errors. 


1.3  Previous  Work 


Two  previous  studies  were  reviewed  which  deal  with  developing  a  transition 
matrix  which  propagates  an  error  vector  to  a  condition  rather  than  through  a  fixed 
time  of  flight.  TTiis  section  briefly  describes  each  one. 

Shepperd  develops  a  method  for  directly  determining  the  six  dimensional 
transition  matrix  corresponding  to  any  conic  state  extrapolation  problem  [6].  His 
method  assumes  a  general  form  for  the  variation  of  the  conic  constraint  equation 
with  three  unknown  functional  coefficients.  He  then  derives  the  state  transition 
matrix  in  terms  of  these  three  unknown  functions,  which  are  different  for  each  conic 
problem.  His  structure  is  formulated  in  terms  of  Goodyear’s  universal  variables. 
After  presenting  a  completely  general  method  for  computing  the  state  transition 
matrix  for  an  arbitrary  conic  extrapolation  problem,  he  demonstrates  the  derivation 
of  the  three  functions  for  the  most  commonly  encountered  single  craft  constraints; 
extrapolation  through  a  fixed  time,  or  through  a  central  angle,  and  extrapolation  to 
a  terminal  radius.  The  advantage  of  this  method  is  that,  in  general,  Kepler’s 
equation  need  not  be  solved,  nor  is  it  necessary  to  determine  the  time  transition 
matrix. 

Tempelman  takes  a  different  approach,  starting  with  the  six  dimensional  time 
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transition  matrix  and  then  allowing  the  time  of  flight  to  vary  [10].  By  assuming 
linear  motion  at  the  terminal  point,  the  vehicle  is  constrained  to  either  a  traverse  or 
a  terminal  cutoff  condition  with  small  variations  in  the  final  time.  An  associated 
time  partial  which  defines  the  change  in  the  time  of  flight  is  developed  for  each 
condition.  Tempelman’s  method  is  also  completely  general,  making  use  of  sensitivity 
vectors,  which  are  the  partial  derivatives  of  the  desired  cutoff  condition  with  respect 
to  the  final  and  initial  state  vectors.  He,  too,  demonstrates  the  derivation  of  the 
sensitivity  vectors  for  several  common  single  craft  traverse  and  terminal  constraints. 
Tempclman’s  approach  is  the  starting  point  for  this  study  because  of  its  inherent 
simplicity  and  versatility  in  handling  a  multitude  of  problems.  Tlie  development  of 
Section  2.3.1  parallels  his  work  in  ref.  [10]. 


1.4  Thesis  Overview 

This  study  develops  a  generalized  condition  transition  matrix  for  a  different 
class  of  constraints;  maneuver  conditions  concerned  with  the  relative  position  of  a 
chaser  and  target  spacecraft  during  rendezvous. 

Chapter  2  presents  the  analytic  development.  First,  the  time  transition  matrix 
is  augmented  in  three  ways: 

•  Time  is  added  to  the  state,  so  that  the  time  partial  is  incorporated  into 
the  final  condition  transition  matrix. 

•  The  state  is  expanded  to  include  both  the  chaser  and  target  vehicles. 

•  Tlie  perturbation  state  is  further  expanded  to  included  navigation  errors 
and  guidance  dispersions. 

The  augmented  condition  transition  matrix  is  then  derived  with  the  help  of  the 
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sensitivity  vector.  Finally,  the  derivation  of  the  sensitivity  vector  is  demonstrated  for 
several  elevation  angle  conditions. 

A  simulation  is  presented  in  Chapter  3,  which  compares  the  analytic  results 
to  a  stochastic  Monte  Carlo  simulation.  The  simulation  was  written  to  validate  the 
analytic  development  and  to  gain  some  insight  into  the  size  of  perturbations  that  can 
be  introduced  without  exceeding  the  region  of  linearity  at  the  terminal  point.  The 
key  components  of  the  simulation  (the  error  generator,  Monte  Carlo,  and  the 
Analytic  routine)  are  described  in  some  detail. 

Chapter  4  presents  the  results  of  the  simulation.  Key  findings  are  made  to 
support  the  validity  of  the  analytic  development  and  the  linearity  assumption.  Two 
coordinate  transformations  are  used  to  show  the  effect  of  the  linearity  assumption 
on  the  final  perturbations  and  to  help  ’see’  what  the  constrained  perturbations  ’look’ 
like  at  the  maneuver  condition. 

Conclusions  drawn  from  the  results  of  the  simulation  are  summarized  in 
Chapter  5.  Lessons  learned  from  this  study  are  presented,  and  topics  for  future 
research  are  mentioned. 
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CHAPTER  2 

ANALYTIC  DEVELOPMENT 


For  certain  rendezvous  mission  architectures,  when  the  primary  goal  is  to 
bring  tw'o  spacecraft  together,  the  actual  time  that  certain  events  or  conditions  occur 
takes  on  less  importance  than  the  fact  that  they  do  occur  at  some  point  in  the 
relative  trajectory.  In  this  case,  the  time  of  flight  of  a  perturbed  trajectoiy  is  allowed 
to  vary  from  the  nominal  to  ensure  that  the  condition  is  still  met  in  the  presence  of 
errors.  While  several  different  conditions  may  be  utilized  during  the  course  of  a 
rendezvous,  the  primary  condition  of  interest  in  this  study  is  the  elevation  angle  of 
the  line  of  sight  to  the  target  with  respect  to  either  the  chaser’s  local  horizontal  or 
its  velocity  vector.  This  provides  a  unique  scalar  condition  to  be  met  at  the  terminal 
condition  point. 

This  chapter  takes  these  conditions  and  develops  an  analytic  method  of 
propagating  a  covariance  matrix  to  satisfy  any  such  condition.  Section  2.1  shows 
how  the  covariance  matrix  is  propagated.  Section  2.2  defines  the  individual  members 
of  the  state,  and  Section  2.3  develops  the  condition  transition  matrix. 
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2  1  Propagating  to  a  Condiaon 


llie  normal  method  of  covariance  matrix  propagation  uses  the  time  transition 
matrix,  <I>  j..  Just  as  this  transition  matrix  takes  an  error  vector  through  a  fixed  time 

of  flight,  it  can  be  used  to  propagate  the  covariance  matrix  through  the  same  fixed 
time  of  flight: 

=  <i>  j  (2.1) 

where  E  is  the  covariance  matrix  and  the  subscripts  /  and  /  designate  the  initial  and 
final  points  on  the  nominal  trajectory. 

The  same  idea  works  for  propagating  to  some  other  scalar  condition.  If  we 
can  find  a  tune  varying  condition  transition  matrix,  which  takes  an  error  vector 

to  a  given  condition  while  allowing  the  time  of  flight  to  vary,  the  same  pre-  and 
post-multiplication  will  propagate  the  initial  covariance  matrix  to  this  condition: 

Deriving  a  method  to  find  this  condition  transition  matrix  is  the  primary  goal  of  this 
chapter. 


2.2  Defining  the  State  and  the  Time  Transition  Matrix 

Before  we  can  find  a  condition  transition  matrix,  we  must,  first,  define  the 
elements  of  the  state  vector  we  wish  to  propagate.  I'his  will  also  lead  to  the 
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appropriate  time  transition  matrix,  another  necessary  ingredient.  The  starting  point 
is  the  traditional  six  dimensional  state  of  one  craft  including  three  components  of 
position  and  three  of  velocity.  The  errors  in  this  state  are  then  defined  as: 
6x  =  [6r  ,  6v]. 


2.2.1  Adding  Time  to  the  State 

Because  we  are  going  to  allow  time  to  slip  to  meet  the  terminal  condition, 
we  need  to  keep  track  of  this  slip.  The  easiest  way  to  do  this  is  to  add  time  to  the 
state.  This  time  ’error’  will  be  defined  as  the  difference  between  the  time  the 
perturbed  state  achieves  the  condition,  and  the  time  the  nominal  state  achieves  the 
condition.  In  most  cases,  the  initial  time  error  will  be  zero,  but  it  is  conceivable  to 
have  an  initial  error  resulting  from  a  prior  time  slip  or  simply  a  small  difference 
between  the  actual  and  nominal  state  update  times. 

For  a  first  look,  the  equation  for  propagating  errors  with  the  six  dimensional 
time  transition  matrix  is  used,  assuming  a  constant  time  of  flight.  The  equation  for 
propagating  the  six  dimensional  state  and  a  simple  equation  for  time  can  then  be 
combined  to  define  a  larger  state  and  time  transition  matrix. 


bXf  = 

bXj 

<1>  0 
^  r  '^6 

bx^ 

^tf  = 

0/  1 

where  0^  is  a  six  dimensional  null  vector.  Notice  that  any  initial  time  error  is  just 

passed  through  to  the  final  state.  In  this  case,  because  the  time  of  flight  is  kept 
constant,  there  is  no  change  in  the  time  error  and  this  error  has  no  effect  on  the  rest 
of  the  state. 
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Taking  this  formulation  one  step  further,  a  small  time  slip,  6tj,  is  introduced 


by  linearizing  the  motion  at  the  terminal  point. 


5x^  =  +  x,6t,* 

bt^  =  5r.  +  6r; 


6x, 

0 

■6x, 

0^  1 

1 

(2.4) 


where  x  is  the  first  derivative  of  the  nominal  state  (referred  to  as  the  dynamical 
state  by  Tempelman).  For  conic  motion,  x  =  [v  ,  a],  the  craft’s  velocity  and 

acceleration  vectors,  with  a  =  -\irlr^.  Tlie  time  slip  is  arbitrary  for  the  moment, 
but  will  provide  the  degree  of  freedom  needed  to  satisfy  some,  as  yet  unspecified, 
scalar  constraint  condition. 

However,  for  rendezvous,  two  state  vectors  must  be  considered. 
Furthermore,  not  only  would  we  like  to  have  the  two  craft  get  together  or  at  least 
arrive  at  a  specified  relative  condition,  but  we  need  this  event  to  occur  at  the  same 
time  for  both  craft.  That  means  the  final  time  error  for  the  chaser  must  equal  the 
final  time  error  for  the  target.  Since  the  two  craft  could  have  different  initial  time 
errors,  we  need  some  way  to  ensure  they  end  up  at  the  condition  with  the  same 
clock  time.  While  we  could  plan  on  having  a  different  time  slip  for  each  craft  to 
account  for  any  initial  time  differences,  an  easier  and  more  straight-forward  solution 
is  to  convert  initial  time  errors  into  initial  position  and  velocity  errors.  Then, 
because  this  initial  time  error  is  ’absorbed’  by  the  rest  of  the  initial  state,  it  makes 
no  contribution  to  the  final  time  error.  As  a  result,  an  equal  time  slip  for  both  craft 
to  meet  the  terminal  condition  requirement  will  give  us  equal  final  time  errors,  as 
required. 

In  ref.  [9],  Tempelman  shows  how  this  might  be  done.  The  basic  idea  is  that 
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since  the  time  slip  is  arbitrary  at  this  point,  no  generality  is  lost  if  a  bias  is 
introduced.  Thus,  defining: 


(2.5) 


where  St,  is  a  new  time  slip  and  Sr,  is  the  bias.  Introducing  this  change  of  variables 
into  equation  (2.4)  yields: 


bXf  =  <Dr5x.  +  if  (bt,  -  6/.) 

btf  =  5^  +  (6^,  -  K 


+ 

- "1 

O 

1 

(2.6) 


This  result  is  basically  the  same  as  equation  (2.4)  except  that  this  initial  time  slip 
error  has  been  absorbed  into  the  state  vector  part  of  the  extrapolation.  Instead  of 
propagating  the  initial  state  errors  a  fixed  time  of  flight  and  leaving  them  at  a 
perturbed  time,  the  errors  are  propagated  by  the  same  time  of  flight  and  then 
’flown’  back  to  the  nominal  time.  Then,  in  both  cases,  they  are  propagated  to 
whatever  time  slip  is  required  to  meet  the  terminal  condition.  Because  the  final 
time  error  is  equal  to  the  time  slip,  both  craft  will  be  on  the  same  clock. 

Note  the  seven  dimensional  time  transition  matrix  is  now  singular,  which  only 
means  we  cannot  go  backward  in  time  and  separate  the  initial  time  error  from  the 
initial  position  and  velocity  errors.  They  are  joined  for  good.  Fortunately,  this  is 
not  a  problem. 


2.2.2  Two  Vehicles 

At  this  point,  we  are  confronted  with  the  task  of  rendezvousing  two 
spacecraft,  each  with  its  own  uncertain  state,  modeled  as  a  nominal  stale  and  a 
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covariance  matrix.  Much  of  current  literature  simplifies  the  problem  by  either 
assuming  perfect  knowledge  of  the  target,  or  grouping  the  uncertainty  of  both  craft 
into  the  chaser’s  covariance  matrix.  The  time  slip  is  then  propagated  along  a 
relative  trajectory  which  is  the  difference  between  the  dynamical  states  of  the  chaser 
and  target  (x^  -  ~  case,  the  effect  is  the  replacement  of  two 

il  state  vectors  by  the  single  relative  state  vector  which  leads  to  suboptimal 
performance. 

In  the  Mars  environment,  with  errors  likely  to  be  larger,  this  approximation 
may  not  be  acceptable.  We,  therefore,  keep  all  the  information  available  by 
augmenting  the  state  to  include  both  vehicles. 


bx,  =  ^  y.  bx.  -  i  +  i  bt. 


e^  Cf  Cf  s 


bx,  =  <5  y  bx,  -  X,  bt,  +  X,  bt, 
7  't  U  7  V  7  * 


bt.  =  bt.  =  6/ 


fix, 

^A  " 

t 

1 

1 

^  T, 

Kt  ^T,  % 


0l  Oj 

aT  ftT 


0 

0 


Oi 

0 

0 


(2.7) 


where  the  subscripts  c  and  t  designate  chaser  and  target,  and  A  designates  the 
augmented  states  and  time  transition  matrix.  Note  that  the  basic  equation  is 
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unchanged. 


2.23  Errors  and  Dispersions 

Covariance  analysis  often  involves  the  use  of  two  basic  types  of  perturbations. 
Navigation  errors,  or  just  errors,  are  the  difference  between  the  navigation,  or 
estimated,  state  and  the  actual  state  of  the  vehicle.  These  errors  develop  from 
imperfect  navigation  measurements  which  are  used  by  a  Kalman  filtering  system  to 
estimate  position  and  velocity.  When  formulated  as  a  covariance  matrix,  they 
represent  the  uncertainty  of  the  navigation  estimate.  Dispersions,  on  the  other 
hand,  are  the  difference  between  the  vehicle’s  actual  state  and  the  nominal,  or 
planned,  state.  They  are  the  result  of  imperfect  maneuver  execution  and  the  use  of 
linearized  guidance  laws  or  imperfect  gravity  models  during  a  mission.  As  the 
vehicle  tries  to  proceed  along  its  nominal  trajectory,  small  errors  in  the  time  of 
execution  of  a  maneuver  or  the  direction  or  duration  of  a  burn  cause  the  vehicle  to 
deviate  from  the  nominal  trajectory.  In  a  covariance  matrix  formulation,  the 
dispersions  represent  the  space  about  a  point  on  the  nominal  trajectory  where  the 
vehicle  should  actually  be.  The  estimated  (navigation)  dispersion  from  nominal  is 
the  sum  of  the  state  error  and  the  actual  dispersion. 


•'mu 

+ 

bx 

= 

bx 

+ 

bx 

= 

+ 

e 

where  x  is  the  six  dimensional  state  of  position  and  velocity  w'ith  subscripts 
designating  the  nominal,  true,  and  navigation  states,  6x  is  now  the  dispersion,  e  is 


Figure  2.1  Perturbed  Position  Vectors 
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the  error,  and  6i  is  the  estimated  dispersion  from  nominal.  Figure  2.1  illustrates 
the  relationship  for  the  position  vector.  Notice  that  dispersions  affect  both  the 
navigation  and  true  state,  while  errors  affect  only  the  navigation  state.  Because 
errors  and  dispersions  are  small  perturbations  about  the  nominal,  they  are  both 
propagated  with  the  same  time  transition  matrix, 

A  time  slip,  or  initial  time  error,  is  converted  to  a  dispersion,  as  developed 
in  Section  2.2.1.  A  navigation  time  error  is  conceivable,  though;  in  ref.  [1],  Battin 
develops  a  procedure  for  tracking  clock  inaccuracies  by  including  time  in  the  state. 
Also,  if  a  spacecraft’s  clock  differed  from  the  ground,  the  time  tag  associated  with 
a  ground  update  would  introduce  further  navigation  errors.  However,  the  accuracy 
of  today’s  clocks  makes  these  errors  small  enough  to  ignore.  This  study  will  assume 
perfect  clocks  and,  therefore,  end  up  with  a  26  dimensional  perturbation  state:  state 
errors,  state  dispersions,  and  time  for  each  of  the  chaser  and  target  vehicles. 
Table  2.1  summarizes  the  elements  of  the  perturbation  state. 

In  constructing  the  final  augmented  system,  it  is  important  to  remember  that 
the  time  slip  affects  dispersions  only,  and  not  errors  (at  least,  to  first  order). 
TTierefore,  the  dynamical  state  components  associated  with  the  navigation  error  in 
the  augmented  dynamical  state  and  time  transition  matrix  are  zero.  Equation  (2.9), 
on  page  33,  shows  the  final  linearized  model  of  conic  motion  for  propagating  the 
errors  and  dispersions  of  two  space  craft  to  a  nominal  final  time  plus  an  equal  time 
slip  for  both  craft.  The  linearized,  variable  time  of  flight  model  of  conic  motion, 
developed  here,  will  be  the  starting  point  for  the  calculation  of  the  condition 
transition  matrix  in  Section  2.3. 
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Table  2.1  Twenty-Six  Dimensional  Perturbation  State  Definition 


Perturbation  Type 

Symbol 

Description 

Navigation  Errors 

Sr, 

Chaser  position  (three  components) 

8v, 

Chaser  velocity  (three  components) 

Sr, 

Target  position  (three  components) 

Sv, 

Target  velocity  (three  components)  | 

Dispersions 

8x, 

Sr, 

Chaser  position  (three  components)  | 

Sv, 

Chaser  velocity  (three  components)  | 

8x, 

Sr, 

Target  position  (three  components)  | 

Sv, 

Target  velocity  (three  components) 

Chaser  time 

St, 

Target  time 
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=  ^T*<: 

e.  = 

7 

fix,  =  4>  r  fix,  -  X  fif  +  X  fir 
fix,  =  ^  fix,  -  X  fir  +  X  fir 

V  7  0  7  * 


fir,  =  fir,  =  fir, 

v  * 


fix,  = 


X 

= 

f 

1 

1 

(2.9) 


<I> 


<I> 


^  */ 


T, 


0 
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2.3  The  Condition  Transition  Matrix 


As  stated  previously,  the  condition  transition  matrix  differs  from  the  time 
transition  matrix  in  that  the  time  of  flight  is  not  held  constant.  Instead,  the  flight 
time  is  allowed  to  vary  so  that  a  scalar  termination  condition  can  be  maintained  in 
the  presence  of  perturbations.  The  methodology  developed  in  this  section  is 
completely  general  and  can  be  used  to  calculate  a  transition  matrix  for  any  scalar 
termination  condition  which  can  be  defined  by  the  state.  This  study  focuses  on  a 
line  of  sight  elevation  angle  which  can  be  used  to  trigger  an  initial  rendezvous  burn. 

Also  discussed  previously,  the  major  simplification  which  allows  this 
development  is  an  assumption  of  linear  motion  at  the  terminal  point.  Because  the 
actual  motion  is  not  linear,  we  are  limited  to  ’small’  time  slips.  The  Monte  Carlo 
simulation,  developed  in  Chapter  3,  will  not  only  validate  the  analytic  development, 
but  also  determine  the  region  of  linearity,  giving  us  a  feel  for  what  is  ’small’. 

2.3.1  Varying  Time  of  Flight  to  meet  a  Scalar  Condition 

Any  arbitraiy  scalar  maneuver  condition,  5,  that  can  be  defined  by  the 
terminal  state,  can  be  used  to  constrain  a  transition  matrix  by  requiring  the  total 
derivative  of  the  condition  to  be  zero. 

6^  =  -^bx.  =  =  0  (2.10) 

dx^  •'  * 

where  k^,  the  .sensitivity,  or  constraint,  vector,  is  the  derivative  of  the  desired 
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condition  with  respect  to  the  terminal  state.  It  should  be  mentioned  that 
Tempelman  allows  for  the  possibility  of  a  constraint  condition  involving  both  initial 
and  final  state  perturbations.  These  are  his,  so  called,  traverse  conditions,  and,  as 
a  result,  his  equations  are  slightly  more  complicated.  Traverse  constraints  can  be 
defined  for  single  vehicle  motion,  but  with  two  vehicle  motion,  as  in  the  case  of 
rendezvous,  they  appear  to  have  little  meaning.  For  this  reason,  this  study  is 
restricted  to  terminal  constraint  conditions,  that  is,  constraints  that  are  functions  of 
the  terminal  state  only. 

'ITie  final  state  is  constrained  by  combining  the  linearized,  variable  time  of 
flight  model  of  coiiic  motion  (from  equation  (2.4)), 

bx^  =  ^  jbx.  +  (2.11) 


with  equation  (2.10),  which  relates  the  desired  constraint  to  the  required  change  in 
flight  time: 

65  =  k'^bxj  =  k'^<b^bx^  +  k’^Xjbt^  =  0  (2-12) 

Notice,  from  equation  (2,12),  the  magnitude  of  k^  is  not  important,  only  its 
direction. 

Equation  (2.12)  can  now  be  solved  for  this  time  slip. 


-k^ibr^x. 
bt  =  - ^ 

*  k^x. 


(2.13) 


and  inserted  into  equation  (2.11): 
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X  k'^ 

6xc  = 

Xy  k 


if  k^ 

I  f 


xjk 


<I>y6x, 


(2.14) 


=  <f>j<l>rbx^ 


- 


which  represents  the  constrained,  linearized  model  of  conic  motion  which  terminates 
on  the  desired  condition.  It  is  completely  general  for  any  state  that  fully  defines  the 
desired  terminal  condition. 

Using  the  augmented  model  developed  in  the  previous  section  yields  a 
desired  advantage.  Because  the  time  error  is  included  in  the  augmented  state,  and, 
f?om  equation  (2.9),  =  6ty,  equation  (2.13)  has  not  been  lost.  It  has  become  a 

part  of  equation  (2.14). 

Equation  (2.14)  shows  that  the  condition  transition  matrix  is  related  to  the 
time  transition  matrbe  by  the  dynamical  state  and  the  sensitivity  vector. 


<Pc  = 


/  - 


Xjk 


xJk 


<b. 


=  A  A 


(2.15) 


Thus,  the  condition  transition  matrix  is  seen  to  be  the  product  of  two  matrices:  an 
idempotent  matrix,  ^  (an  idempotent  matrix  is  any  matrix  A,  such  that  A^=A),  and 
the  time  transition  matrix  [10]. 
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The  result  in  equation  (2.14)  can  be  rewritten  to  show  that  the  final 
constrained  perturbation  is  the  product  of  the  idempotent  matrix  and  final 
perturbation  at  the  nominal  terminal  time. 


!>Xc  =  ^ 


(2.16) 


Because  this  matrix  is  idempotent,  it  can  be  interpreted  as  a  ’shaping’  matrix, 
shaping  the  perturbations  to  allow  the  final  state  to  conform  to  the  required 
condition.  The  idempotent  property  confirms  that  once  a  vector  is  shaped,  any 
attempt  to  shape  it  again  with  the  same  matrix  will  have  no  effect;  it  is  already 
constrained  to  the  proper  condition. 

Combining  equation  (2.15)  with  equation  (2.2)  and  then  equation  (2.1),  yields 
the  desired  propagation  of  the  covariance  matrix  to  a  condition  and  shows  how  the 
idempotent  matrix  ’shapes’  a  covariance  matrbe  which  has  been  propagated  to  a 
nominal  final  time. 


Ec  =  9  c  E,  <1^1 

=  (2-17) 

2.3.2  Maneuver  Conditions  and  the  Sensitivity  Vector 

While  the  methodology  developed  is  completely  general  and  can  be  used  to 
calculate  a  transition  matrix  for  any  terminal  maneuver  condition,  this  study  focuses 
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on  the  elevation  angle  of  the  line  of  sight  from  the  chaser  craft  to  the  target  craft. 
This  elevation  angle  is  an  observable  condition  which  might  be  used  to  trigger  a 
rendezvous  burn.  Four  definitions  for  the  elevation  angle  are  developed.  The  first 
is  the  elevation  of  the  line  of  sight  with  respect  to  the  local  horizontal  of  the  chaser 
craft,  and  the  second  is  the  same  angle  with  a  correction  if  the  target  and  chaser  are 
not  co-planer.  The  last  two  elevation  angles  are  definea  relative  to  the  chaser’s 
velocity  vector.  Figure  2.2  illustrates  the  elevation  angles  with  no  out  of  plane 
correction. 

When  the  elevation  angle  is  defined  with  respect  to  the  local  horizontal  of 
the  chaser  craft,  it  can  be  written  as  a  function  of  the  position  vectors  of  the  two 
spacecraft. 


sin4>j 


\^r\\rj 


(2.18) 


where  Ar  =  r,  -  is  the  relative  position  vector.  The  corre*.  ion  made  in  0, 

removes  the  out  of  plane  component  of  the  relative  position  vector,  essentially 
projecting  the  target’s  position  onto  the  orbital  plane  of  the  chaser  craft. 


sin<J), 


I'l'-  - 

Ar^  r. 


(2.19) 


|Ar  -  -L(r>,)AJ|rJ 


where  h  =  r  x  v  is  the  orbital  angular  momentum  vector.  4>2  is  also,  therefore,  a 
function  of  the  chaser’s  vehx'iiy  vector  unless  the  tw'o  craft  are  co-planar,  in  which 
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Figure  2.2  Flevation  Angle  Conditions 


case  equation  (2.19)  simplifies  to  equation  (2.18). 

(f)j  and  (f>^  are  defined  similarly,  but  with  respect  to  the  chaser’s  velocity 

vector: 


cos  4)3 


(2.20) 


and 


sin  4)^ 


Ar  -  —(r^h  )h 

T 

V, 

.t  \  *  tj  e 

C 

Ar^ 


'i'-  - 


(2.21) 


The  sensitivity  vector,  which  is  the  derivative  of  the  desired  condition  with 
respect  to  the  final  state,  is  a  row  vector  containing  the  partial  derivatives  of  the 
condition  with  respect  to  each  member  of  the  state: 


d4>  d4)  d4>  d4>  d4>  34> 
de^  ’  de,  ’  dx^  *  dx/  dt/  dr, 


(2.22) 


where 


de  dx 


dr*  dv 


(2.23) 


Any  members  of  the  state  which  are  not  a  direct  function  of  the  condition  have  a 
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zero  partial  derivative; 


=  0  ^  _  0T 

dt^  dt,  dv,~^ 


(2.24) 


for  all  elevation  angles  defined. 

Deriving  the  remaining  partials  for  each  elevation  angle  is  a  straight-forward, 
but  messy  exercise  in  differentiation.  A  few  important  steps  and  the  results  of  the 
differentiation  are  contained  in  Appendix  A.  Note  the  magnitude  of  the  sensitivity 
vector  is  not  needed  in  calculating  the  condition  transition  matrix,  so  any  scalar 
values  common  to  all  partials  can  be  neglected,  if  desired. 
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CHAPTER  3 
SIMULATION 


A  Monte  Carlo  simulation  was  written  to  validate  the  analytic  development 
of  Chapter  2  and  to  gain  some  knowledge  of  the  size  of  perturbations  that  can  be 
introduced  without  exceeding  the  region  of  linearity  at  the  terminal  point.  This 
chapter  describes  the  implementation  of  this  simulation.  Section  3.1  gives  an 
overview  of  the  main  components  of  the  simulation  and  describes  two  ’off  the  shelf 
procedures  used.  Section  3.2  shows  how  representative  error  vectors  were  generated 
from  a  covariance  matrbc,  detailing  the  statistical  theory  involved.  The  central 
procedure  that  propagated  all  the  perturbed  states  to  the  condition  and  calculated 
the  statistics  is  the  subject  of  Section  3.3,  and  finally,  the  implementation  of  the 
analytic  solution  is  discussed  in  Section  3.4. 


3.1  Simulation  Overview 

The  Monte  Carlo  simulation  has  six  main  components,  as  detailed  in 
Table  3.1  (Tlieir  hierarchy  is  shown  in  Figure  3.1).  Simulation  Setup,  the  initializer 
and  main  executive,  and  two  ’off  the  shelf  procedures,  Telev  and  Kepler,  are 
addressed  in  this  section.  The  remaining  components;  the  error  generator,  Monte 
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Figure  3.1  Hierarchy  of  Simulation  Components 
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Table  3.1  Major  Components  of  Computer  Simulation 


Program  Component 

Purpose 

Simulation  Setup 

Initializes  simulation  and  acts  as  main 
executive. 

Elevation  Angle  Search 

Finds  the  time  a  desired  elevation  angle 
condition  will  occur. 

Kepler 

Propagates  a  spacecraft’s  state  to  a  new  time 
and  computes  the  time  transition  matrix. 

Error  Generator 

Generates  a  statistically  representative  set  of 
error  vectors  from  an  initial  covariance  matrix. 

Monte  Carlo 

Propagates  the  perturbed  states  to  the 
condition  and  calculate^  the  st^tisliVs 

Analytic 

Implements  the  analytic  results  of  Chapter  2. 

Carlo,  and  Analytic;  are  the  subject  of  the  following  sections. 

Inputs  to  the  simulation  are: 

•  Initial  nominal  state  of  the  chaser  in  inertial  frame. 

•  Initial  nominal  state  of  the  target  in  inertial  frame. 

•  Initial  nominal  time. 

•  Desired  elevation  angle. 

•  First  guess  for  nominal  time  to  achieve  elevation  angle  condition. 

•  Type  of  elevation  angle  (from  the  four  definitions). 

•  Initial  augmented  covariance  matrix  containing  error  and  dispersion 
covariance  data  for  the  chaser  and  target. 

•  Covariance  matrix  scale  factor. 

Outputs  are: 

•  Final  twenty-six  dimensional  inertial  covariance  matrix  propagated  to 
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the  desired  condition  by  the  Monte  Carlo  routine. 

•  Final  twenty-six  dimensional  inertial  covariance  matrix  propagated  to 
the  desired  condition  by  the  Analytic  routine. 

•  Three  dimensional  Monte  Carlo  estimated  position  dispersion 
covariance  matrix  rotated  to  a  relative  line  of  sight  frame. 

•  Three  dimensional  Analytic  estimated  position  dispersion  covariance 
matrix  rotated  to  a  relative  line  of  sight  frame. 

•  Three  dimensional  Monte  Carlo  true  position  dispersion  covariance 
matrix  rotated  to  a  relative  line  of  sight  frame. 

•  Three  dimensional  Analytic  true  position  dispersion  covariance  matrix 
rotated  to  a  relative  line  of  sight  frame. 

One  sigma  errors  are  also  calculated  from  the  main  diagonals  of  each  matrix. 

Simulation  Setup  initializes  the  simulation  by  solving  the  nominal  rendezv'ous 
problem.  It  calls  Telev  to  find  the  nominal  final  time  for  the  desired  elevation  angle 
condition,  and  Kepler  to  propagate  the  chaser  and  target  states  to  this  nominal 
condition  (see  Figure  3.2).  It  then  acts  as  the  executor,  calling  the  error  generator 
to  obtain  the  initial  augmented  error  set  and  then  Monte  Carlo  and  Analytic  to 
generate  the  final  condition  covariance  matrices. 


^  to  condition 

^  X  '  t 

nom^  Telev,  Kepler 

nom^  ’  nom 
c  c 

Figure  3.2  Propagating  the  Nominal  State  to  an  Elevation  Angle  Condition 


Kepler  is  an  ’off  the  shelf  procedure  written  by  Sian  Shepperd  which 
propagates  a  position  and  velocity  vector  and  computes  a  six  dimensional  time 
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transition  matrix  for  a  given  time  of  flight.  The  procedure  solves  Kepler’s  problem 
and  calculates  the  transition  matrix  using  Goodyear’s  universal  variables.  Further 
information  on  the  Kepler  procedure  can  be  found  in  ref.  [7]. 

The  elevation  angle  search  component  is  also  an  ’off  the  shelf  procedure 
called  Telev  and  taken  from  the  space  shuttle’s  on  orbit  software.  It  uses  an 
iterative  routine  to  find  the  time  a  desired  elevation  angle  condition  will  occur  and 
has  been  adapted  to  use  the  elevation  angles  defined  in  eqs.  (2.18)  through  (2.21). 
It  was  set  to  use  tvs'o  body  orbital  mechanics.  Further  information  on  Telev  can  be 
found  in  ref.  [5]. 

3.2  Error  Generator 

A  procedure  was  written  to  generate  any  number  of  error  vectors  which 
accurately  represent  a  given  twenty-four  dimensional  covariance  matrix  (Initial  time 
errors  were  assumed  to  be  zero,  with  both  craft  starting  on  the  same  clock).  The 
error  generator  takes  advantage  of  the  ’square  root’  covariance  matrix  devised  by 
Jim  Potter,  as  cited  by  Battin  [1]. 

The  development  takes  three  steps  [2]; 

•  A  vector,  y,  of  independent,  gaussian  random  numbers  with  unit 
variance  has  covariance, 

E[yyT]  =  ^  =  / 
where  /  is  the  iflentity  matrix. 

•  TTie  nxn  symmetric,  positive  semidefinite  covariance  matrix,  P, 
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representing  the  expected  value  of  the  outer  product  of  an  n 
dimensional  error  vector,  e,  with  itself, 


P  =  E[ee'^  = 


(3.2) 


can  be  written  as  the  product  of  a  matrix,  W,  and  its  transpose: 

p  =  (3.3) 


W  is,  therefore,  a  square  root  of  P. 

•  A  vector, 

X  =  Wy  (3-4) 

will  have  covariance, 

E[xjcT  =  W{E\yy'^)W'^  =  =  P 

The  square  root  matrix,  W,  is  not  unique.  One  method  of  finding  a  square 
root  matrix  is  to  use  the  eigenvectors  and  eigenvalues: 

p  =  VAV^  =  VA*A»V^  =  WW^  ^  ^ 


where  V  is  the  matrix  of  unit  eigenvectors  (also  called  the  modal  matrix)  and  A  is 
the  diagonal  eigenvalue  matrix. 

One  commonly  used  error  generating  method  takes  advantage  of  this  fact 
by  interpreting  the  eigenvectors  as  ’directions’  in  rz-space  along  which  the 
components  of  the  error  vectors  are  uncorrelated.  The  eigenvalues  are  the  variances 
of  these  components.  The  desired  error  vectors  are  then  generated  using  the 
following  equation: 
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(3.7) 


i«l 

where  r  is  an  independent,  gaussian  random  number  with  unit  variance,  and  A,,  and 
V,  are  the  i’'  eigenvalue  and  eigenvector  of  P  [4].  Equation  (3.7)  can  be  seen  to  be 
equivalent  to: 

X  =  =  Wy 


Calculating  the  eigenvectors  and  eigenvalues  of  a  twenty-four  dimensional  matrix  is 
a  time  consuming  process,  though.  Fortunately,  its  also  avoidable. 

Because  any  square  root  matrix  will  do  the  trick,  a  much  easier  method 
assumes  a  lower  triangular  form  for  W  [1].  This  assumption  allows  the  desired 
square  root  matrix  to  be  determined  by  the  straightforward  solution  of  a  series  of 
simultaneous  algebraic  equations.  While  there  are  a  variety  of  ways  of  doing  this, 
the  following  easily  programmed,  recursive  algorithm  for  an  nxn  positive 
semidefinite  matrix,  P,  was  borrowed  from  Cholesky  and  Banachiewicz,  as  cited  in 
ref.  [Ij: 

For  /=!,  2,  ...,  n,  the  elements  of  W  are  calculated  from 


w,.  = 


Pa  - 


y-i 


0 


Pji  - 


for  j<  1 

for  y  =  i  +  l,j  +  2,...,n 


(3.9) 


where  w  and  p  are  the  elements  of  W  and  P. 
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Because  of  the  efficiency  of  this  technique,  it  was  also  used  to  improve  the 
computer’s  gaussian  pseudo-random  number  generator  when  generating  ’small’ 
vector  sets.  After  generating  a  set  of  independent,  pseudo-random  vectors,  z,  their 
covariance  was  calculated  to  be  close  to  /,  but  not  exact  (the  diagonal  elements 
ranged  from  0.92  to  1.07,  and  the  off  diagonal  elements  were  as  large  as  0.08  for 
1,001  twenty-four  dimensional  vectors).  The  computed  covariance, 

Efzz'n  =  zz^  =  B=  UU^  -  7 

where  B  is  the  near-identity  covariance  matrix,  and  U  is  its  square  root,  was, 
therefore,  ’divided’  out  to  produce  a  set  of  vectors  that  really  had  independent 
elements  with  unit  variance: 


y  =  U^z 

E[y>1  =  C/'Ve[zz'T){/''^  =  U'BU^  =  / 


(3.11) 


The  new  random  vectors,  y,  indeed,  had  a  covariance  equal  to  the  identity  (to  at 
least  15  figures).  And  the  computed  covariance  of  the  resulting  error  vector  set, 
X  =  Wy,  matched  the  desired  covariance,  P,  just  as  well.  This  ’fix’  was  derived 
independently,  but  was  later  found  to  be  previously  developed  by  Suddath  in  ref. 
[8]- 


3.3  Monte  Carlo 


The  Monte  Carlo  component  was  written  to  solve  the  rendezvous  problem 
for  many  randomly  perturbed  states  obtained  from  the  initial  augmented  covariance 
matrix.  The  Monte  Carlo  component  propagates  the  generated  error  vectors  to  the 
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final  condition,  then  re-calculates  the  statistics,  forming  the  final  condition 
covariance  matrix. 

For  each  augmented  error  vector,  Monte  Carlo  builds  an  initial,  true  and 
navigation  state  from  the  initial  nominal  state,  as  defined  in  Chapter  2,  for  the 
chaser  and  the  target: 


(3.12) 


For  each  navigation  state,  it  calls  Telev  to  find  the  perturbed  time  for  the  desired 
elevation  angle  condition,  and  Kepler  to  propagate  the  estimated  chaser  and  target 
states  to  this  perturbed  condition.  It  then  calls  Kepler  again  to  propagate  the  true 
chaser  and  target  states  to  the  same  perturbed  time  (see  Figure  3.3). 


X 

to  condition 

X  ■  t 

nav . 

% 

Telev,  Kepler 

>  nav^ 

X. 

tmjie . 

to  t 

nav  Q 

-►  X, 

Kepler 

true^ 

% 

c 

Figure  3.3  Propagating  the  Perturbed  State  to  an  Elevation  Angle  Condition 

The  final  chaser  and  target  condition  perturbations  arc  then  calculated  from  these 
states  by  re-writing  equation  (3.12): 
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=  t 

-  t 

c 

nofHc 

naVf- 

bXf. 

-  X 

-  X  ^ 

nofe 

MMf. 

bXr, 

= 

-  X 

c 

tnu^ 

Komf. 

«c  = 


where  the  subscript,  C,  designates  a  state  or  perturbation  propagated  to  the  desired 
elevation  angle  condition. 

Finally,  an  augmented  final  condition  perturbation  vector,  is  built 

as  defined  in  Table  2  1.  These  augmented  vectors  are  then  used  to  calculate  the 
final  condition  covariance  matrix: 


(3.14) 


3.4  Analytic  Method 

If  the  previous  section  seemed  complex,  it  should  serve  to  evoke  an 
appreciation  for  an  analytic  method  of  propagating  a  covariance  matrix  to  a 
maneuver  condition.  TTie  simulation’s  Analytic  component  simply  computed  one 
(twenty-six  dimensional)  condition  transition  matrix  using  the  equations  developed 
in  Chapter  2  and  summarized  here  for  clarity.  As  throughout  this  study,  the  six 
dimensional  time  transition  matrix  is  assumed  to  be  known,  and  was  generated  by 
Kepler. 
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The  first,  and,  perhaps,  most  difficult,  step  is  to  calculate  the  augmented 
sensitivity  vector, 


where 


and 


T  _  8<t) 

_  d(t>  d4>  d(|>  d<J>  d4>  84> 

'dT/  aT’  aT’  ^ 

Cf  tf  tf  t,  c  t 


(3.15) 


a^)  _  a4>  _  a<j)  a4> 

de  dx  dr'  dv 


(3.16) 


^  =  0 
ar,  ■  at,  ■ 


av, 


(3.17) 


The  remaining  partials  are  located  in  Appendix  A. 

The  augmented  dynamical  state  and  time  transition  matrix  are  then 
computed  as  follows: 


where 


and 


(3.18) 
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<I> 


T, 


-X. 


T. 


-X. 


0 


(3.19) 


which  provides  all  the  pieces  necessary  to  compute  the  condition  transition  matrix: 


<I>. 


J^f 


•T. 


(3.20) 


Individual  augmented  perturbation  vectors  were  then  propagated  for 
comparison  to  Monte  Carlo, 

bXf,  =  (3.21) 


and  the  final  condition  covariance  matrix  was  computed  in  one  step: 


Ec= 


(3.22) 


The  format  for  the  output  covariance  matrices  is  presented  in  Chapter  4. 
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CHAPTER  4 
RESULTS 


Several  test  cases  were  used  to  validate  the  results  of  Chapter  2,  investigate 
the  region  of  linearity  at  the  terminal  point,  and  get  a  feel  for  the  geometry  of  the 
perturbations  at  the  desired  condition.  The  results  of  these  test  cases  are  presented 
in  this  chapter. 

The  problem  setup,  including  the  initial  covariance  matrix  and  the  nominal, 
initial  and  final  vehicle  states  are  presented  in  Section  4.1.  Section  4.2  compares  the 
final  condition  covariance  matrices  generated  by  the  Monte  Carlo  and  Analytic 
components  of  the  simulation,  verifying  they  agree.  In  Section  4.3,  these  matrices 
are  transformed  into  a  curvilinear  coordinate  frame  to  bring  out  the  effect  of  the 
linearity  assumption.  Finally,  in  Section  4.4,  an  effort  is  made  to  visualize  the 
perturbations  at  a  maneuver  condition  by,  first,  using  a  different  coordinate 
transformation  to  ’look’  at  the  results  in  a  relative,  line  of  sight  frame  to  see  if  the 
errors  collapse  along  the  line  of  sight,  and  then  by  comparing  the  estimated 
elevation  angle  error  to  the  true  elevation  angle  dispersion. 
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direction 
of  motion^ 


geocentric  radius 
vector  of  vehicle 


vehicle 
orbital  plane 

NAME:  U,  V,  W  coordinate  system  [5J. 

ORIGIN:  Point  of  interest. 

ORIENTATION:  The  U-V  plane  is  the  instantaneous  orbit  plane  at  epoch.  The  U 
axis  lies  along  the  geocentric  radius  vector  to  the  vehicle  and  is 
positive  radially  outward. 

The  W  axis  h'es  along  the  instantaneous  orbital  angular 
momentum  vector  at  epoch,  and  is  positive  in  the  direction  of  the 
angular  momentum  vector. 

The  V  axis  completes  the  right-handed  orthogonal  system. 

CHARACniRIS  nCS:  Quasi-inertial,  right-handed,  Cartesian  coordinate  system.  This 

system  is  quasi-incrtial  in  the  sense  that  it  is  treated  as  an  inertial 
coordinate  system,  but  it  is  redefined  at  each  point  of  interest. 


Rgure  4.1  UVW  Coordinate  System 
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4.1  Problem  Setup 


All  test  cases  used  the  same  initial  covariance  matrix.  This  covariance 
matrix  reflects  the  uncertainty  that  exists  in  the  ground  uplink  of  a  shuttle  and  target 
state  during  a  shuttle  rendezvous  mission.  The  full  six  by  six  chaser  and  target 
covariance  matrices,  listed  in  Table  4.1,  are  expressed  in  the  UVW  coordinate  frame 
(see  Figure  4.1),  and  were  rotated  to  the  inertial  frame  for  the  initialization  of  each 
test  case.  Four  assumptions  were  made  in  forming  the  initial  augmented  covariance 
matrix: 

•  Uncertainty  in  the  errors  and  dispersions  of  each  craft  are  equal. 

•  Correlations  between  chaser  and  target  craft  are  zero. 

•  Correlations  between  errors  and  dispersions  are  zero. 

•  Initial  time  errors  are  zero. 

The  augmented  covariance  niatrbc  was  then  multiplied  by  a  scale  factor  to  analyze 
the  effects  of  scaled  perturbations  on  the  results. 

Two  different  rendezvous  orbits  were  used.  The  first  had  the  target  craft  in 
a  200  NM  altitude  circular,  equatorial  orbit  with  the  chaser  ten  miles  below  in  a 
co-circular,  co-planar  orbit.  The  nominal  final  condition  point  was  planned  to  occur 
with  the  chaser  craft  on  the  positive,  inertial  X  axis.  The  second  case  had  the  target 
in  a  200  by  29400  mile  elliptic,  equatorial  orbit  (eccentricity  =  0.8)  with  the  chaser, 
again,  ten  miles  below  in  a  co-elliptic,  co-planar  orbit.  Periapsis  coincided  with  the 
positive,  inertial  X  axis,  and  the  nominal  final  condition  point  was  planned  to  occur 
with  the  chaser  craft  on  the  positive  Y  axis  (true  anomaly  =  90°).  The  initial 
nominal  state,  as  input  to  the  simulation,  and  the  final  nominal  state,  as  propagated 
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Table  4.1  Initial  Covariance  Matrices  in  UVW  Frame 


Chaser’s  Error  and  Dispersion  Covariance  Matrix: 


3.83193E+05 

-5.60173E+04 

1.34268E+03 

4.66259E+02 

-3.02118E+02 

-7.65116E-01 


-5.60173E+04 

2.35247E+06 

-4.99017E+03 

-2.18982E+03 

-4.19257E+01 

9.47873E-01 


1.34268fi+03 

-4.99017E+03 

1.17616E+06 

6.96687E+00 

-2.i5600E+00 

-4.06157E+02 


4.66259E+02 

-2.18982E+03 

6.96687E+00 

2.57924E+00 

-2.96924E-01 

-1.98502E-03 


-3,021 18E+ 02 
-4.19257E+01 
-2.15600E+00 
-2.96924E-01 
2.47009E-01 
1.22858E-03 


-7.651 16E-01 
9.47873E-01 
-4.06157E+02 
-1.98502E-03 
1.22858E-03 
3.81924E-01 


Target’s  Error  and  Dispersion  Covariance  Matrix: 


1.35124E+04 

-1.48492E+04 

3.59579E+02 

2.47387E+01 

-1.55673E+01 

7.76503E-02 


-1.48492E+04 
9.08787E+05 
- 1. 47445  E +04 
-1.07590E+03 
1.32557E+01 
-5.09445E+00 


3.59579E+02 
- 1. 47445  E+ 04 
1.06319E+06 
1.86755E+01 
-2.08800E+00 
5.60670E+02 


2.47387E+01 
-1.07590E+03 
1.86755E+01 
1.28142E+00 
-2.35343E-02 
6.0494  IE -03 


-1.55673E+01 

1.32557E+01 

-2,08800E+00 

-2.35343E-02 

1.82250E-02 

-9.91980E-04 


7.76503E-02 
-5.09445E+00 
5.60670E+02 
6.0494  IE- 03 
-9.91980E-04 
4.46224E-01 
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to  the  requested  elevation  angle  by  Simulation  Setup,  are  listed  in  Table  4.2.  The 
format  used  in  this  table  to  print  the  fourteen  dimensional  state  vector  is  as 
follows:^ 


r.  1 


The  desired  elevation  angle  for  each  case  was  varied  according  to  the  orbit 
and  the  elevation  type  requested,  so  the  nominal  final  points  for  all  elevation  types 
would  be  the  same.^  For  the  co-circular  orbits,  a  desired  angle  of  25.08  degrees  was 
used.  For  the  co-elliptic  orbits,  when  the  elevation  was  measured  relative  to  the 
chaser’s  local  horizontal  (elevation  types  1  and  2),  51.82°  was  used;  when  the 
elevation  was  measured  relative  to  the  chaser’s  velocity  vector  (types  3  and  4),  13.15° 
was  used. 


4.2  VaHdating  the  Analytic  Result 


The  analytic  result  of  Chapter  2  was  validated  in  three  ways.  The  first  way 
used  a  separate  simulation  which  propagated  small  deterministic  dispersions  with  a 
fourteen  dimensional  condition  transition  matrix  (no  navigation  errors).  The  Kepler 
procedure  was  then  used  to  propagate  the  initial  perturbed  state  to  the  same  final 
time  determined  by  the  analytic  method.  The  final  elevation  angle  of  the  perturbed 


’The  units  of  measure  used  for  all  results  are;  feet,  feet  per  second,  seconds,  and  degrees,  as 
appropriate. 

’Reference  to  elevation  type  in  the  results  is  from  the  four  elevation  angle  definitions  of  Section  23.2. 
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Table  4.2  Nominal  Simulation  States 


Initial  Nominal  State  for  Co-Circular  Case: 
-3.84059n+06  -2.17811E+07  0.00000 

-3.56721  E-t-06  -2.18891E+07  0.00000 

-1.53011E+03  - 1.5301  lE-l-03 

2.48447E-I-04 

2.48654E+04 

-4.38080E-(-03 

-4.05225E-^03 

0.00000 

0.00000 

Final  Nominal  State  for  Co-Circular  Case: 
2.21171E-»-07  2.95695E-08  0.00000 
2.21774E+07  1.29025E-f-05  0.00000 

0.00000  0.00000 

-4.72937E-11 

-1.46569E+02 

2.52280E-t-04 

2.51930E-I-04 

0.00000 

0.00000 

Initial  Nominal  State  for  Co-Elliptic  Case: 
-1.83697E+08  -3.40423E-(-07  0.00000 
-1.83700E-t-08  -3.43293E-H07  0.00000 
- 2.3 1422E -1-04  -2.31422E-I-04 

3.42607E-t-03 

3.44948E-H03 

-  3. 44055  E -1-03 
-3.43606E-I-03 

0.00000 

0.00000 

Final  Nominal  State  for  Co-Elliptic  Case: 
-1.93715E-07  3.98168E-t-07  0.00000 

-2.17846E-I-05  4.00938E+07  0.00000 

0.00000  0.00000 

-1.88024E-I-04 

-1.87778E-I-04 

1.50471  E+04 
1.49204E-I-04 

0.00000 

0.00000 
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states,  as  produced  by  the  analytic  method  and  Kepler,  were  compared  to  the 
desired  elevation  angle.  The  results  of  this  comparison,  for  several  simple 
dispersions,  are  presented  in  Table  4.3.  The  analytic  method  is,  indeed,  finding  the 
correct  time  slip  and  is  propagating  the  dispersions  to  the  desired  condition. 

The  remaining  methods  of  validation  used  the  simulation  described  in 
Chapter  3.  TTie  first  of  these  was  a  direct  comparison  of  the  augmented  condition 
covariance  matrices  produced  by  the  Analytic  and  Monte  Carlo  components.  In  the 
interest  of  conciseness,  only  the  ’one  sigma’  uncertainties  (the  square  root  of  the 
main  diagonals)  are  presented  here.  For  comparison,  the  uncertainties  of  the  initial 
covariance  matrix  for  the  co-circular  case  are  shown  in  Table  4.4.  The  format  used 
in  all  remaining  tables  to  print  the  twenty-six  one  sigma  uncertainties  is  as  follows: 


] 

The  order  is  as  defined  in  Table  2.1,  with  a  position  and  velocity  perturbation  vector 
on  each  row. 

Table  4.5  compares  the  results  for  a  co-circular  case  in  which  the  initial 
covariance  matrix  was  scaled  by  a  factor  of  one  sixteenth.  As  hoped,  both 
simulation  components  are  producing,  essentially,  the  same  condition  covariance 
matrices.  Also  noteworthy,  is  the  relative  sizes  of  the  perturbations.  Tlie  errors 
stayed  about  the  same  size  as  the  initial  errors,  while  the  dispersions  grew' 
dramatically.  This  is  due  to  the  time  slip.  The  average  time  slip  (about  17  seconds, 
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Table  4.3  Elevation  Angle  Accuracy  Using  Analytic  Method 


Initial  Error 

St,  (s) 

Elevation  Error  (deg) 

Analytic 

Kepler 

5r,  =  10  ft 

0.193 

-1.2E-6 

-2.4E-7 

Sr,  =  10  ft 

0.804 

-2.0E-5 

-8.9E-6 

SVjj  =  0.1  ft/s 

4.573 

-9.2E-4 

-5.9E-4 

Svy  =  0.1  ft/s 

-4.602 

-5.9E-4 

-2.3E-4 

Sr.  =  [1,  -3,  0] 

Sv.  =  [-0.1,  0.3,  0] 

Sr,  =  [-1,  4,  0] 

Sv,  =  [-0.04,  -0.09,  0] 

St.  =  lE-8 

St,  =  0 

-7.691 

-2.3E-3 

-1.4E-3 

Table  4.4  Initial  Covariance  Matrix  -  Co-Circular,  Covariance  Scale  Factor  1,  Inertial 
Frame 


One  Sigma  Uncertainties  From  Main  Diagonal: 


1.52061E+03 

9.43578F+02 

1.52061E+03 

9.43578E+02 

0.00000 


6.50703E+02 

1.78773E+02 

6.50703E+02 

1.78773E+02 

0.00000 


1.08451E+03 
1.031  llE+03 
1.084.S1E+03 
1.03111E+03 


6.472 16E-01 
2.41615E-01 
6.472 16E -01 
2.41615E-01 


1.55157E+00 

1.11412E+00 

1.55157E-t-00 

1.11412E-I-00 


6.18000E-01 

6.68000E-01 

6.18000E-01 

6.680001:-01 
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Table  4.5  Comparison  of  Final  Condition  Covariance  Matrices  -  Co-Circular, 
Elevation  Type  1,  Covariance  Scale  Factor  0.0625,  Inertial  Frame 


One  Sigma  Uncertainties  From  Main  Diagonal: 

2.39691E+02 

2.41656E+01 

2.39691E-I-02 

2.55629E+03 

1.74177E+01 

7.69378E+02 
2.8641  lE+02 
4.38899E+05 
4.38855E+05 
1.74177E+01 

ANALYTIC: 

1.66196E+02  8.00467E-01 

1.12794E+02  3.12385E-01 

1.66196E+02  5.00695E+02 

1.12794E+02  4.98525E+02 

2.25876E-01 

2.24063E-02 

2.25876E-01 

2.89725E+00 

2.89095E-01 

3.11802E-01 

2.89095E-01 

3.11802E-01 

2.41163E+02 

2.47189E+01 

7.31356E+03 

7.95173E+03 

1.73170E+01 

7.69301  E+02 
2.8635  lE+02 
4.36274E+05 
4.3623  lE+05 
1.73170E+01 

MONTE  CARLO: 

1.66159E+02  8.00686E-01 

1.12647E+02  3.12309E-01 

1.66004E+02  4.97700E+02 

1.12990E+02  4.95544E+02 

2.25955E-01 
2.33777E-02 
8.37151  E+00 
9.03020E+00 

2.89123E-01 

3.11871E-01 

2.89239E-01 

3.11710E-01 

Table  4.6  Comparison  of  Final  Condition  Covariance  Matrices  - 
Elevation  Type  3,  Covariance  Scale  Factor  0.01,  Inertial  Frame 

Co-Elliptic, 

One  Sigma  Uncertainties  From  Main  Diagonal: 

4.91245E+03 

3.6.S232E+03 

1.13589E+06 

1.13830E-I-06 

6.05487E+01 

4.28278E+03 

3.15071E+03 

9.09298E+05 

9.04527E+05 

6.05487E+0I 

ANALYTIC: 

5.92013E+02  3.65586E-01 

6.72483E+02  1.49546E-01 

5.92013E+02  3.65586E-01 

6.72483E+02  2.91747E+00 

2.28709E+00 

1.68834E+00 

5.36496E+02 

5.30792E+02 

2.82168E-01 

3.00228E-01 

2.82168E-01 

3.00228E-01 

4.90978E+03 
3.65n.80E+03 
1.13306E+06 
1.1 35481-: +06 
6.042411-: +01 

4.28722E+03 
3.15846E+03 
9.06152E+05 
9.01415E+05 
6.0424  lE+01 

MONTE  CARLO: 
5.92390E+02  3.72.361E-01 

6.71716E+02  1.54584E-01 

5.92372E+02  1.34829E+01 

6.73257E+02  1.34324E+01 

2.291 96E  + 00 
1.69562E+00 

5  34592E+02 
5.28922E+02 

2.82249E-01 

3.00885E-01 

2.82243E-01 

3.00131E-01 
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in  this  case)  shows  itself  in  the  dispersions,  while  the  errors  are  unaffected,  as  stated 
in  the  development. 

The  co-elliptic  case,  presented  in  Table  4.6,  validates  the  development  using 
the  elevation  angle  as  measured  from  the  chaser’s  velocity  vector.  Again,  Analytic 
and  Monte  Carlo  produce,  essentially,  the  same  results.  In  this  case,  the  elevation 
type  makes  a  significant  difference.  As  shown,  with  elevation  type  3,  the  average 
time  slip  is  60.5  seconds.  When  the  same  case,  with  the  same  covariance  scale 
factor,  was  run  with  the  elevation  defined  off  the  local  horizontal  (elevation  type  1, 
not  shown),  the  average  time  slip  required  was  only  24.1  seconds  (As  mentioned  in 
the  previous  section,  the  requested  elevation  angle  was  varied  for  each  elevation 
type  so  the  nominal  final  positions  of  both  craft  would  be  the  same). 

The  final  method  of  validation  looked  at  the  average  elevation  error.  The 
sensitivity  vector  and  the  final  condition  transition  matrbc  are  used  to  calculate  a 
scalar  variation  in  the  elevation  angle: 


5<t> 


34) 


"  *7*  A 


(4.1) 


Of  course,  this  value  is  guaranteed  to  be  (and  is)  zero  using  the  analytic  method, 
because 


(4-2) 

but  the  fact  it  is  near  zero  when  using  Monte  Carlo’s  output  validates  both  the 
sensitivity  vector  derivation  and  the  Monte  Carlo  method.  Table  4.7  summarizes  the 
average  elevation  error  (the  square  root  of  equation  (4.1))  for  several  test  cases. 
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Table  4.7  Estimated  Average  Elevation  Angle  Error  Calculated  From  tlie  Condition 
Covariance  Matrix 


Orbit 

Elevation 

Type 

Covariance 
Scale  Factor 

Average  Elevation  Error 

Analytic 

Monte  Carlo 

co-circular 

1 

25 

2E-5 

2.9 

4 

7E-7 

0.6 

1 

2E-6 

0.2 

6.25E-2 

lE-6 

lE-2 

co-elliptic 

1 

lE-2 

2E-7 

lE-2 

3 

lE-2 

3E-6 

4E-2 

4.3  The  Linearity  Assumption 

The  assumption  of  linear  motion  at  the  terminal  point  cannot  be  taken  for 
granted,  because,  obviously,  it  must  break  down  somewhere.  For  ’small’  enough 
perturbations  in  time,  though,  a  linear  region  must  exist.  To  bring  out  the  effect  of 
the  linearity  assumption,  the  individual  twenty-six  dimensional  perturbation  vectors, 
which  were  propagated  to  the  condition  by  Monte  Carlo  and  Analytic,  were 
transformed  to  a  local  vertical,  curvilinear  (LVC)  coordinate  system  (see  Figure  4.2). 
These  transformed  vectors  were  then  used  to  compute  a  new  LVC  condition 
covariance  matrix.  By  separating  the  downrange  and  altitude  (Z,,,) 

perturbations,  a  vehicle,  which  has  been  propagated  linearly,  can  be  seen  ’gaining 
altitude’  above  the  actual  curved  flight  path.  The  linearity  assumption  is,  therefore. 
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Figure  4.2  Nominal  Vehicle-Centered,  Local  Vertical,  Curvilinear  Coordinate  System 
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seen  to  introduce  an  altitude  dispersion  that  is  not  really  there.  For  small  time  slips, 
this  added  dispersion  is  also  small,  but  for  larger  slips,  the  dispersion  gets  to  be 
excessive. 

The  co-circular  orbit  case,  with  two  different  sized  initial  covariance  matrices 
is  presented  to  show  this  effect.  In  Table  4.8,  the  same  case  and  covariance  scale 
factor  as  for  Table  4.5  are  used.  Note  the  altitude  dispersions  for  the  analytic 
component  are  larger  than  for  Monte  Carlo.  For  Monte  Carlo,  the  full  effect  of  the 
time  slip  is  seen  in  the  downrange  dispersions.  For  Table  4.9,  the  covariance  scale 
factor  is  increased  to  1,  and  the  effect  is  more  pronounced.  The  altitude  dispersion 
has  grown  to  about  one  tenth  the  downrange  dispersion  in  the  analytic  covariance 
matrix. 

Next,  with  a  good  knowledge  of  the  direct  effect  of  the  linearity  assumption, 
it  is  necessary  to  measure  the  accuracy  of  the  analytic  component  in  keeping  the 
elevation  error  small  as  the  size  of  the  initial  perturbations  are  increased.  A  more 
direct  measure  of  the  actual  estimated  average  elevation  error  was  obtained  by 
calculating  the  elevation  angle  of  the  final  navigation  states  produced  by  the 
Analytic  and  Monte  Carlo  components  from  each  initial  perturbation  vector  and 
then  calculating  how  this  varied  from  the  nominal.  Table  4.10  presents  this  result 
for  the  same  cases  presented  in  Table  4.7.  Note  that  when  the  elevation  error  is 
calculated  by  this  method,  Monte  Carlo  is  very  consistent  (reflecting  the  elevation 
tolerance  required  of  the  iterative  elevation  angle  search  procedure),  but  Analytic 
starts  to  break  down  with  very  large  perturbations.  The  effect  of  large  perturbations 
on  the  ability  of  the  sensitivity  vector  to  calculate  an  accurate  elevation  angle 
variance  can  now  be  recognized  in  Table  4.7,  too. 

Another  measure  of  the  accuracy  of  the  analytic  method  is  its  ability  to  find 
the  correct  time  slip.  Table  4.11  shows  the  one  sigma  time  slips  obtained  by 
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Table  4.8  Comparison  of  Final  Condition  Covariance  Matrices  -  Co-Circular, 
Elevation  Type  1,  Covariance  Scale  Factor  0.0625,  LVC  Frame 


One  Sigma  Uncertainties  From  Main  Diagonal: 


7.69021E-H02 

2.86382E-I-02 

4.38733E-I-05 

4.38699E+05 

1.74177E-t-01 


1.66196E-»-02 

1.12794E-r02 

1.66196E-r02 

1.12794E-(-02 

1.74177E+01 


ANALYTIC: 


2.39948E-r02 

2.38315E+01 

7J8400E+03 

7J4485E+03 


4.97599E-01 

4.96460E-02 

4.97534E-01 

4.96366E-02 


2.89095E-01 

3.n802E-01 

2.89095E-01 

3.11802E-01 


1.13777E-01 

3.70283E-02 

1.14011E-01 

3.70405E-02 


7.69801  E +02 
2.86489E+02 
4.36355E+05 
4.36324E+05 
1.73170E+01 


1.66159E+02 

1.12647E+02 

1.66004E+02 

1.12990E+02 

1.73170E+01 


MONTE  CARLO: 
2.39560E+02  4.97720E-01 

2.30692E+01  4.96883E-02 

2J9518E+02  4.97626E-01 

2.30616E+01  4.96737E-02 


2.89123E-01 

3.11871E-01 

2.89239E-01 

3.11710E-01 


1.I2898E-01 

3.70330E-02 

1.12962E-01 

3.70350E-02 


Table  4.9  Comparison  of  Final  Condition  Covariance  Matrices  -  Co-Circular, 
Elevation  Type  1,  Covariance  Scale  Factor  1,  LVC  Frame 


One  Sigma  Uncertainties  From  Main  Diagonal: 


ANALYTIC; 

9.99848E+02 


1.-30122E+02 

1.17176E+05 

1.I6779E+05 


1.97034E+00 

1.98017E-01 

I.96838E+00 

1.97890E-01 


3.04724E+03 
1.13879E+03 
1.74521E+06 
1.74514E+06 
6.967  lOE +01 


3.08807E+03 
1.14568E+03 
1.75.33  IE +06 
1.75318E+06 
6.95810E+01 


6.64784E+02 

4.51177E+02 

6.64784E+02 

4.51177E+02 

6.96710E+01 


6.66729E+02 
4.5.3618E+02 
6.64571  E+02 
4. 5  9074  E  + 02 
6.95810E+01 


MONTE 
9.54406E+02 
9.27440E+01 
9.53523  E+02 
9.27005E+01 


CARLO: 
1.98230E+00 
1.998.52i:-01 
1.980.331-4  00 
1. 99803  li-Ol 


1.15638E+00 

1.24721E+00 

1.15638E+00 

1.24721E+00 


1.154Q3E  +  00 
1.24606E  +  00 
1. 1 56.55  E  + 00 
1.24348E+00 


4.98496E-01 

1.48148E-01 

5.0.3737E-01 

1.48286E-01 


4.57677E-01 
1.47775E-01 
4.59001  E-01 
1.47745E-01 
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Table  4.10  Estimated  Average  Elevation  Angle  Error  Calculated  From  the  Error 
Vectors 


Orbit 

Elevation 

Covariance 

Average  Elevation  Error 

Type 

Scale  Factor 

Analytic 

Monte  Carlo 

co-circular 

1 

25 

2.7 

4E-3 

4 

0.6 

4E-3 

1 

0.2 

3E-3 

6.25E-2 

lE-2 

4E-3 

co-elliptic 

1 

lE-2 

2E-2 

5E-3 

3 

lE-2 

4E-2 

3E-3 

Table  4.11  Average  Time  Slip  Comparison 


Orbit 

Elevation 

Type 

Covariance 

Scale 

Factor 

Average  Time  Slip 

Average 
Time  Slip 
Difference 

Analytic 

Monte  Carlo 

co-circular 

1 

25 

348 

335 

13 

4 

139.3 

138.7 

0.7 

1 

69.7 

69.6 

0.1 

6.25E-2 

17.4 

17.3 

0.1 

co-elliptic 

1 

lE-2 

24.1 

24.3 

-0.2 

3 

lE-2 

60.5 

60.4 

0.1 
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Analytic  and  Monte  Carlo,  and  their  difference.  The  average  time  slip  obtained 
from  Analytic  starts  to  diverge  from  that  of  Monte  Carlo  as  the  perturbations  grow. 
Both  measures,  the  elevation  angle  error  and  the  average  time  slip  difference,  give 
a  feel  for  what  ’small’  perturbations  are  for  these  orbits. 

Unfortunately,  the  definition  of  ’small’  will  change  for  each  rendezvous 
scenario.  Factors  which  affect  how  quickly  altitude  errors  build  up,  and,  thus,  how 
large  an  acceptable  time  slip  can  be,  are  the  curvature  of  the  actual  flight  path  and 
the  speed  of  the  vehicle.  This  means  the  size  of  the  orbit  and,  for  elliptic  orbits,  the 
vehicle’s  position  in  the  orbit  will  affect  the  definition  of  small.  Many  factors  will 
also  affect  how  quickly  the  condition  constraint  is  changing,  and,  therefore,  the  size 
of  the  time  slip  required.  In  the  final  analysis,  though,  the  size  of  the  time  slips 
required  is  primarily  a  function  of  how  well  the  problem  is  defined.  It  may  well  be 
the  case  that  the  analytic  method  developed  here  will  not  be  applicable  to  certain 
ill-defined  problems  requiring  large  time  slips  which  violate  the  region  of  linearity. 


4.4  Visualizing  Perturbations  at  a  Maneuver  Condition 

Trying  to  think  in  twenty-six  dimensions  is  a  little  difficult.  This  makes 
’seeing’  what  the  perturbations  at  a  maneuver  condition  look  like  equally  difficult. 
Any  effort  to  reduce  the  dimensionality  of  the  problem  is  a  great  help.  The  LVC 
frame  was  useful  to  show  that  the  perturbations  were  really  about  the  same  size,  but 
just  spread  up  and  down  range  from  nominal.  Another  transformation  to  a  relative, 
line  of  sight  (LOS)  coordinate  system  is  also  helpful  (see  Figure  4.3  for  a  definition). 
By  going  to  a  relative  frame,  which  expresses  the  chaser’s  position  relative  to  the 
target,  the  dimensionality  is  cut  to  twelve  (there  is  no  relative  time  difference 
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NAME:  Line  of  sight  coordinate  system. 

ORIGIN:  Chaser  vehicle  center  of  mass. 


ORIENTA'IION:  The  X|o,-Y,o,  plane  is  the  instantaneous  plane  spanned  by  the 

geocentric  radius  vectors  of  the  chaser  and  target  vehicles.  The 
X|a,  axis  lies  along  the  relative  position  vector  and  is  positive 
toward  the  target  vehicle.  This  is  the  line  of  sight  from  the  chaser 
to  the  target. 

The  Z,o,  axis  is  normal  to  the  X|o,-Y,o,  plane,  and  is  positive  in  the 
direction  of  a  vector  which  is  the  vector  cross  product  of  the 
chaser  and  target  geocentric  radius  vectors. 

The  Y^,  axis  is  perpendicular  to  the  relative  position  vector  and 
completes  the  right-handed  orthogonal  system. 

CIIARAC'll' RISTICS:  Right-handed,  Cartesian,  rotating  coordinate  system.  Tliis  system 

is  designed  so  the  elevation  angle,  as  defined  in  equation  (2.18), 
will  lie  in  the  X,o,-Y,o,  plane.  A  slightly  different  line  of  sight 
frame  is  used  with  each  elevation  angle  definition  to  keep  the 
angle  in  this  plane. 


Figure  4.3  Line  of  Sight  Coordinate  System 
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between  vehicles).  Neglecting  velocity  errors  and  dispersions,  we  can  then  look  at 
estimated  position  dispersions  and  true  position  dispersions  separately,  and  the 
problem  is  down  to  three  dimensions.  The  former  shows  how  the  navigation 
position  estimate  is  dispersed  from  the  nominal,  and  the  latter  shows  the  true 
position  dispersions  from  the  nominal. 

Intuition  would  suggest  that  when  the  line  of  sight  is  held  constant  relative 
to  the  local  horizontal  or  the  velocity  vector,  that  the  error  perpendicular  to  the  line 
of  sight  and  in  the  relative  plane  (Y,<^)  would  be  zero  when  the  condition  is 
achieved.  The  error  vector  is  said  to  ’collapse’  along  the  line  of  sight.  This  is 
impossible  to  see  in  an  inertial  frame,  though,  because  the  local  horizontal  and  the 
velocity  vector,  and  thus  the  LOS  frame,  are  rotating  with  time.  And  because  time 
is  allowed  to  slip,  the  times  associated  with  the  nominal  and  perturbed  states  are 
different.  Figure  4.4  illustrates  how  perturbing  a  state  from  the  nominal  with  a  time 
slip  affects  the  LOS  frame,  and,  therefore,  shows  that  when  time  is  allowed  to  slip, 
the  LOS  frame  for  the  perturbed  state  must  also  be  rotated  to  compensate. 

An  inertial,  time  varying  perturbation  vector  or  covariance  matrix  cannot, 
then,  be  simply  rotated  to  the  LOS  frame.  In  fact,  a  different  perturbation  vector 
is  needed.  In  an  inertial  frame, 

but  to  obtain  a  time  varying,  LOS  perturbation  vector,  the  LOS  frame  for  the 
perturbed  state  must,  first,  be  rotated  to  account  for  the  time  slip: 

where  is  a  rotation  matrix, 
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A(bibt^)  = 


(4.5) 


cos(wfifp  sin(o>fifp  0 
-sin((o6r^)  cos((o6t^)  0 
0  0  1. 

and  0)  is  the  angular  velocity  of  the  chaser  orbit.  This  is  the  only  rotation  needed 
when  the  line  of  sight  is  held  constant  relative  to  the  local  horizontal.  When  it  is 
held  constant  relative  to  the  velocity  vector,  an  additional  rotation  must  be  made  to 
account  for  the  angular  rotation  of  the  velocity  vector  with  respect  to  the  local 
horizontal.  This  is  accomplished  by  substituting  (o  +  y)  for  o  in  equation  (4.5), 
where  y  is  the  time  rate  of  change  of  the  flight  path  angle: 

y  =  it  ^  (4.6) 

dt 

where  p  is  the  gravitational  parameter.  These  new  LOS  perturbation  vectors,  which 
were  calculated  using  the  twenty-sbc  dimensional  inertial  perturbation  vectors 
propagated  to  the  condition  by  Analytic  and  Monte  Carlo,  were  then  used  to 
compute  a  new  LOS  condition  covariance  matrix. 

Often,  in  covariance  analysis,  the  question  is,  what  happens  to  the  true  state 
when  maneuvers  are  made  based  on  the  navigation  estimate  of  the  state.  In  this 
case,  the  navigation  state  is  constrained  to  a  condition,  but  the  true  state  is  not.  The 
transformation  to  the  LOS  coordinate  system  lets  us  see  the  effect.  In  Table  4.12, 
the  co-circular  orbits  case  with  a  1/16  covariance  scale  factor  is  presented  in  the 
LOS  frame.  The  estimated  position  dispersion  and  true  position  dispersion 
covariance  matrices,  as  calculated  from  the  inertial  perturbation  vectors  produced 
by  Analytic  and  Monte  Carlo,  are  shown,  as  well  as  their  associated  one  sigma 
uncertainties.  As  suggested  by  intuition,  the  estimated  position  dispersion 
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Table  4.12  Comparison  of  Final  Condition  Covariance  Matrices  -  Co-Circular, 
Elevation  Type  1,  Covariance  Scale  Factor  0.0625,  LOS  Frame 

RELATIVE  ANALYTIC  ESTIMATED  POSITION  DISPERSION  COVARUNCE  MATRIX: 
6.261913020919n+05  5.3473293385 15n-h03  -2.540297561 973 E -1-02 

5.3473293385 15E-(-03  7.03191 185801  lE-t-02  -6.229824767212E-t-01 

-2.540297561973E-f-02  -6.229824767212E-f-01  8.068720629973E-)-04 

One  Sigma  Uncertainties  From  Main  Diagonal: 

7.913225019497E-(-02  2.651775227656E-I-01  2.840549353553E-(-02 

REIATIVE  MONTE  CARLO  ESTINUTED  POSITION  DISPERSION  COVARUNCE  NUTRIX: 
6.138037833267E+05  1.444583188955E+03  3.4593 18429782E+02 

1.444583 188955E-t-03  8.276183539467E-H01  -  3.99928698781 6E +01 

3.4593 18429782E+02  -  3.999286987816E+01  8.062224383743E+04 

One  Sigma  Uncertainties  From  Main  Diagonal: 

7.834563059461  E+02  9.097353208196E+00  2.839405639169E+02 

RELATIVE  ANALYTIC  TRUE  POSITION  DISPERSION  COVARUNCE  NUTRIX: 

1.644869483409E+06  6.355364742566E+05  -2.929294 122991  E+02 

6.355364742566E+05  3.020465862788E+05  -  2.364784855042E+02 

-2.929294 122991  E+02  -2.364784855042E+02  4.034360314986E+04 

One  Sigma  Uncertainties  From  Main  Diagonal: 

1.282524652164E+03  5.49587651 1338E+02  2.008571710193E+02 


RELATIVE  MONTE  CARLO  TRUE  POSITION  DISPERSION  COVARUNCE  MATRIX: 


1.626745360030E+06 

6.2767902621 17E+05 

4.211591 1.50865E+02 

6.2767902621 17E+05 

2.979182713376E+05 

2.452058269794E+02 

4.21 1591 150865E+02 

2.452058269794E+02 

4.036887999681E  +  04 

One  Sigma  Uncertainties  From  Main  Diagonal: 

1.275439281201E:+03 

5.458188997622E+02 

2.00920()836074i:+02 

(navigation’s  deviation  from  nominal)  perpendicular  to  the  line  of  sight  and  in  the 
relative  plane  (Y^)  is  very  small  for  both  Analytic  and  Monte  Carlo,  but  the  true 
position  dispersion  (true  deviation  from  nominal),  in  the  same  direction,  is  not. 
Table  4.13  reveals  the  same  effect  for  the  co-elliptic  case  with  a  0.01  covariance 
scale  factor.  It  really  is  easier  to  ’see’  what  is  happening  with  the  problem  reduced 
to  three  dimensions. 

The  effect  on  the  true  state  caused  by  performing  maneuvers  based  on  the 
estimated  state  can  also  be  seen  by  comparing  the  average  elevation  angle  error 
calculated  from  the  estimated  state  with  the  average  elevation  angle  dispersion 
calculated  from  the  true  state.  The  true  average  elevation  dispersion  can  be 
calculated  in  two  ways,  just  as  the  estimated  average  error  was.  Table  4.7  showed 
the  estimated  average  elevation  error  as  calculated  with  the  sensitivity  vector  and  the 
condition  covariance  matrbc  using  equation  (4.1).  The  true  average  elevation 
dispersion  can  also  be  calculated  using  this  equation  by  substituting  a  sensitivity 
vector  which  reflects  only  the  sensitivity  of  the  elevation  angle  to  dispersions.  This 
is  accomplished  by  nulling  the  first  twelve  components  of  k/.  The  resulting  true 
average  elevation  dispersions  are  presented  in  Table  4.14. 

The  second  method  of  calculating  the  true  average  elevation  dispersion  uses 
the  variance  from  nominal  of  the  true  elevation  angle  as  calculated  from  the  final 
true  states  produced  by  the  Analytic  and  Monte  Carlo  components.  This  is  similar 
to  the  method  used  to  find  the  estimated  average  elevation  errors  in  Table  4  10. 
This  second  method  is  not  subject  to  the  linearity  problems  that  the  sensitivity  vector 
encounters  for  large  perturbations.  TTie  true  average  elevation  errors  calculated  by 
this  method  are  shown  in  Table  4.15. 
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Table  4.13  Comparison  of  Final  Condition  Covariance  Matrices  -  Co-Elliptic, 
Elevation  Type  3,  Covariance  Scale  Factor  0.01,  LOS  Frame 


RELATIVE  ANALYTIC  ESTIMATED  POSITION  DISPERSION  COVARUNCE  MATRIX: 
4.282563658181E+07  -9.308798138141E+04  2.023131415154E-K)4 
-9.308798138141E+04  1.696932409177E+04  -6.360480 197909E+ 03 
2.023 131415 154E-J-04  -6.360480197909E-r03  1.605425422832E-I-06 
One  Sigma  Uncertainties  From  Main  Diagonal: 

6.544 129933139E+03  1  J02663582502E-r02  1.267053835806E-^03 

RELATIVE  MONTE  CARLO  ESTIMATED  POSITION  DISPERSION  COVARUNCE  NUTRIX: 
4.291 740970 190E-t- 07  -9.38949435/ I88E-r04  2.863330763954E-1-04 
-9.389494357 188E-f04  5.0051 1 29641 3 9E 3- 04  -8.108242958814E-^03 
2.863330763954E+04  -8.108242958814E+03  1.601 982277946E+ 06 
One  Sigma  Uncertainties  From  Main  Diagonal: 

6.551 1380463 17E-(-03  2.2372 10978906E+ 02  1.265694385682E+03 

RELATIVE  ANALYTIC  TRUE  POSITION  DISPERSION  COVARUNCE  MATRIX: 
4.552474790435E-t-07  8.876379130440E-^06  9.790509630079E-(-03 
8.876379130440E-I-06  3.32000721 1203E-r06  5.104268670316E+03 
9.790509630079E+03  5.1042686703I6E+03  8.0271271 14158E -1-05 
One  Sigma  Uncertainties  From  Main  Diagonal: 

6.74720296896  lE-f  03  1.822088694659E-1-03  8.959423594271E-I-02 

RELATIVE  MONTE  CARLO  TOTAL  POSITION  DISPERSION  COVARUNCE  MATRIX: 
4.540100903208E-(-07  8.814150151238E-1-06  4.9232 1 231 85 13E-t- 03 
8.814150151238E-(-06  3.3385 15525613E-1-06  -9.8390762171 19E-I-03 
4.923212318513E-t-03  -9.8390762171 19E-1-03  8.05 1403937792E -1-05 
One  Sigma  Uncertainties  From  Main  Diagonal: 

6.7380270875 14E-(- 03  1.827 160508990E-r  03  8.972961572297E-t-02 
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Table  4.14  True  Average  Elevation  Angle  Dispersion  Calculated  From  the 
Condition  Covariance  Matrix 


Orbit 

Elevation 

Type 

Covariance 
Scale  Factor 

Average  Elevation  Error 

Analytic 

Monte  Carlo 

co-circular 

1 

25 

4.4 

6.8 

4 

1.8 

2.1 

1 

0.9 

0.9 

6.25E-2 

0.2 

0.2 

co-elliptic 

1 

lE-2 

0.3 

0.3 

3 

lE-2 

0.3 

0.3 

Table  4.15  True  Average  Elevation  Angle  Dispersion  Calculated  From  the  Error 
Vectors 


Orbit 

Elevation 

Covariance 

Average  Elevation  Error 

Type 

Scale  Factor 

Analytic 

Monte  Carlo 

co-circular 

1 

25 

5.2 

5.1 

4 

1.8 

1.8 

1 

0.9 

0.9 

6.25E-2 

0.2 

0.2 

co-elliptic 

1 

lE-2 

0.3 

0.3 

3 

lE-2 

0.3 

0.3 
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CHAPTER  5 
CONCLUSION 


Future  missions  to  Mars  will  stretch  current  theories  and  practices  in  mission 
planning  and  operations  well  beyond  those  used  to  go  to  the  moon.  The  goal  of  this 
study  has  been  to  develop  a  method  for  propagating  a  covariance  matrix  to  a 
maneuver  condition  to  be  used  in  linear  covariance  analysis  for  planning  the 
rendezvous  phase  of  these  space  missions.  With  the  generalized  formulation  of  a 
condition  transition  matrix,  an  analytic  method  of  propagating  an  augmented 
covariance  matrix  to  any  scalar,  terminal  maneuver  condition  has  been  presented. 
This  augmented  covariance  matrbe  includes  navigation  errors,  dispersions,  and  time 
errors  for  both  the  chaser  and  target  craft.  First,  the  method  was  analytically 
developed,  then  validated  by  comparison  to  a  stochastic  Monte  Carlo  simulation  for 
the  case  of  several  elevation  angle  conditions  which  might  be  used  to  trigger  an 
initial  rendezvous  burn.  The  key  assumption  of  linearity  at  the  terminal  point  was 
substantiated  with  the  same  simulation.  This  chapter  will  summarize  the  key  results 
of  the  simulation  in  Section  5.1,  and  suggest  some  future,  related  areas  for  research 
in  Section  5.2. 
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5.1  Summary  of  Results 


The  primary  result  of  this  study  is  a  valid  analytic  method  of  propagating  to 
a  condition.  Three  important  conclusions  can  be  drawn  from  the  results  of  the 
comparison  of  this  analytic  method  with  a  stochastic  Monte  Carlo  simulation: 

•  The  condition  transition  matrix,  developed  in  Chapter  2,  is  a  valid 
analytic  means  of  propagating  a  covariance  matrix  to  a  maneuver 
condition. 

•  The  assumption  of  linearity  at  the  terminal  point  is  valid  for  ’small’ 
initial  perturbations. 

•  For  an  elevation  aiigle  condition,  the  final  estimated  state  errors  at  the 
condition  do  collapse  along  the  line  of  sight  between  the  spacecraft. 

This  section  will  expand  on  these  three  conclusions. 

The  analytic  '•esult  (the  propagation  of  an  initial  augmented  covariance 
matrix  to  a  maneuver  condition  using  a  variable  timp  of  flight  condition  transition 
matrix)  was  validated,  first,  by  accurately  propagating  deterministic  error  vectors  to 
the  desired  condition.  The  elevation  angle  was  held  constant,  and  the  final  state  was 
verified  with  a  separate  Kepler  propagation  to  the  same  perturbed  time.  The 
analytic  method  was  then  shown  to  produce  the  same  final  condition  covariance 
matrix  as  that  of  the  stochastic  Monte  Carlo  simulation.  In  particular,  the  result  was 
shown  for  a  co-circular  orbital  case  A'ith  the  elevation  angle  referenced  to  the  local 
horizontal,  and  for  a  co-elliptic  case  with  the  elevation  angle  referenced  to  the 
chaser’s  velocity  vector.  Finally,  the  variation  in  the  elevation  angle  was  calculated 
from  the  final  condition  co.'ariance  matrices.  This  variation  was  shown  to  be  small. 
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confirming  that  the  final  state  had  been  constrained  to  the  desired  condition. 

The  linearity  assumption  was  validated  for  ’small’  initial  perturbations.  The 
definition  of  small  was  investigated,  and  as  stated,  will  be  unique  to  the  orbits  of  the 
rendezvousing  spacecraft.  Using  the  co-circular  case  and  a  transformation  to  a  local 
vertical,  curvilinear  coordinate  system,  the  effect  of  the  linearity  assumption  was 
shown  to  be  an  added  altitude  error  in  the  final  dispersions  of  both  craft.  Then,  by 
scaling  the  initial  covariance  mrtrix,  the  linearity  assumption  was  tested  for  various 
sized  errors.  The  analytic  method  proved  to  accurately  calculate  the  necessary  time 
slip  and  achieve  the  desired  elevation  angle  for  an  initial  covariance  matrix 
representing  the  uncertainty  of  a  ground  uplinked  shuttle  and  target  state.  When 
these  uncertainties  were  doubled  (with  a  covariance  scale  factor  of  4),  the  time  slip 
still  agreed  with  the  Monte  Carlo  simulation  to  within  0.7  seconds,  and  the  average 
elevation  angle  error  was  held  to  0.6  degrees.  Only  for  much  larger  perturbations 
(covariance  scale  factor  25),  did  the  linearity  assumption  start  to  break  down. 

A  different  tr^msformation,  to  a  rotating,  line  of  sight  coordinate  system,  was 
necessary  to  show  that,  when  the  elevation  angle  is  constrained,  the  total 
perturbation  perpendicular  to  the  line  of  sight  between  the  craft  will  be  small.  It 
was  shown  that,  oecause  the  LOS  frame  rotates  with  time,  and  time  is  allowed  to 
vary,  a  different  perturbation  vector  was  needed  to  account  for  the  rotated  LOS 
frame  of  the  pertnrbed  state  relative  to  the  LOS  frame  of  the  nominal  state.  Once 
this  LOS  perturbation  or  was  calculated,  its  resulting  covariance  matrix  did 
reveal  that  the  perturbations  had  collapsed  along  the  line  of  sight. 

The  condition  transition  matrix,  developed  in  Chapter  2,  will  analytically 
propagate  a  covariance  matrix  to  a  maneuver  condition.  This  method  will  prove  to 
be  an  effective,  time  .  .ving  tool  for  the  covariance  analysis  necessary  to  the  planning 
of  the  rendezvous  phase  of  any  future  space  mission. 
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5.2  Future  Research 


As  is  always  the  case,  there  is  more  work  yet  to  be  done  than  has  already 
been  accomplished.  Three  areas  of  future  research,  directly  related  to  the  condition 
transition  matrix,  suggest  themselves. 

The  definition  of  what  is  small,  as  it  relates  to  the  linearity  assumption,  was 
investigated  for  two  very’  general  rendezvous  orbit  profiles,  one  co-circular  and  the 
other  co-elliptic.  Because  the  size,  shape,  position,  and  relative  size  and  orientation 
of  the  rendezvous  orbits  all  influenee  the  altitude  error  and  affeet  how  fast  the 
elevation  angle  is  ehanging,  and,  therefore,  how  large  the  perturbations  can  be 
before  the  linearity  assumption  is  invalid,  more  work  is  needed  to  define  what  is 
small  for  specific  rendezvous  scenarios. 

If  it,  then,  proves  that  the  size  of  the  anticipated  errors  and  the  planned 
rendezvous  profile  will  necessitate  time  slips  beyond  the  region  of  linearity  at  the 
terminal  point,  work  might  be  done  to  expand  this  region  of  linearity.  One 
possibility  is  to  include  second  order  effects  in  the  condition  transition  matrix. 
Another  possibility  would  be  to  go  to  a  classical  element  formulation  (or  a  more 
modern  non-singular  element  formulation)  which  avoids  these  linearity  problems. 
Unfortunately,  element  formulations  are  difficult  to  deal  with,  both  algorithmically 
and  computationally.  However,  for  some  problems,  it  may  be  necessary. 

The  final  direct  area  of  research  invc!  es  the  transformations  to  the  LVC  and 
LOS  coordinate  systems.  Both  transformations  are  non-linear,  and  were,  therefore, 
performed  on  the  individual  perturbation  vectors  after  they  were  propagated  by 
Analytic  or  Monte  Carlo.  Covariance  matrices  were  then  calculated  from  these 
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transformed  vectors.  An  Analytic  method  for  directly  transforming  the  inertial 
condition  covariance  matrix  should  be  possible  and  would  further  increase  the 
flexibility  and  value  of  the  condition  transition  matrix. 

Not  mentioned  yet,  are  the  areas  which  go  beyond,  and  were  listed  as 
motivation  for,  this  study.  Armed  with  an  analytic  method  of  propagating  a 
covariance  matrix  to  a  maneuver  condition,  questions  related  to  the  Mars  highly 
elliptic  rendezvous  problem  can  be  more  easily  addressed.  Some  of  these  questions 
are:  the  size  and  shape  of  the  catch-up  orbit,  where  in  the  orbit  the  rendezvous 
should  occur  and  where  it  should  be  initiated,  what  condition  to  trigger  the  intercept 
burn  on,  the  geometry  of  the  final  approach  to  the  target,  and  how  to  make  the 
rendezvous  tolerant  to  navigation  uncertainties.  There  is,  indeed,  more  work  yet  to 
be  done. 
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APPENDIX  A 


This  Appendix  contains  the  results,  and  a  few  helpful  steps,  of  differentiating 
the  four  elevation  angle  definitions  of  Chapter  2  to  obtain  the  elements  of  k^,  the 
sensitivity  vector. 


A.l  Elevation  Type  1 


For  the  angle  between  the  line  of  sight  and  the  local  horizontal, 

Ar'f  r 


lArl  |r 


(A.1) 


two  good  first  steps  are: 


Ar  =  r,  -  r. 


8Ar  _  _  dAr 

dTc  dr, 


and 


2Ar 


Ar^ 

dAr 

dr. 


Ar^Ar 

-2Ar^ 


dAr  dAr 


dr. 


dr. 


All 

lArl 


(A.2) 
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and  the  result: 


(A.3) 


(A.4) 


34), 


a,r. 


|Ar^rJ  Ar^ 


(A.S) 


where 


sec4>j 
Ar’  rl 


(A.6) 
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A.2  Elevation  Type  2 


When  correcting  for  non-co-planar  orbits, 


a  few  more  helpful  steps  are  needed  [1]: 


(A.7) 


e  T  (A.8) 


(A.9) 


(A.  10) 
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Similar  results  are  found  in  ref.  [1]  for  the  partials  with  respect  to  the  velocity 
vector. 

With  a  few  more  steps,  the  desired  partials  are: 
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A.3  Elevation  Type  3 
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A.4  Elevation  Type  4 


leads  to 


sin<t»4 


Ar^v 


K 

Ar’^v 


(A.20) 


«4''c 


(Ar^vJ 


^  Ar^ 

*  (Ar^,)(r;vJ](r7*,f 

^(A>'’'v,)(r7*,)(r,xv.)T 


-  Ar'' 


(A.21) 
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where 


dv. 


=  a. 


-^(r7*,f(Ar''»,)(r;v,)r; 


+  V, 


Ar^ 


*  (A^^v.)U(*=  -  r^3(r7».f 

"t 


-  4r’ 


(A.22) 


5<t>4 

dr. 


=  «4Vc 


(A.23) 


OL,  = 


-CSC  <1)4 


Ar’  -  Mr^h,f 


(A.24) 
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